Large scale behavior and statistical equilibria in rotating flows 
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We examine long-time properties of the ideal dynamics of three-dimensional flows, in the presence 
or not of an imposed solid-body rotation and with or without hehcity (velocity- vorticity correlation). 
In all cases the results agree with the isotropic predictions stemming from statistical mechanics. No 
accumulation of excitation occurs in the large scales, even though in the dissipative rotating case 
anisotropy and accumulation, in the form of an inverse cascade of energy, are known to occur. 
We attribute this latter discrepancy to the linearity of the term responsible for the emergence of 
inertial waves. At intermediate times, inertial energy spectra emerge that differ somewhat from 
classical wave-turbulence expectations, and with a trace of large-scale excitation that goes away for 
long times. These results are discussed in the context of partial two-dimensionalization of the flow 
undergoing strong rotation as advocated by several authors. 



I. INTRODUCTION 

Geophysical and astrophysical flows are turbulent. 
The dimensionless parameter that measures such a state, 
the Reynolds number, is huge for the atmosphere and 
the ocean, as well as for the sun and stars, the interplan- 
etary medium, and the interstellar medium. Nonlinear 
effects, e.g., through advection, have thus ample time to 
act; they develop strong localized intermittent structures 
such as vortex filaments (or current sheets when dealing 
with magnetic fields), and a myriad of modes are excited. 
One thus resorts to a statistical description of such flows 
and examine the behavior of probability density func- 
tions (PDFs) as well as scaling exponents of structure 
functions. At second order, one obtains a power-law for 
the distribution of energy among the Fourier modes, in 
a direct cascade of energy to small scales where the en- 
ergy is finally dissipated [l[. Another striking feature of 
turbulent flows is that, in some cases (including in three 
space dimensions in the presence of rotation) , large-scale 
modes are excited as well 0-0|- 

A deep understanding of turbulence is still lacking, in 
part because of the large number of interacting modes 
due to the convolution in Fourier space stemming from 
the nonlinearities, quadratic in the incompressible case. 
Dimensional analysis and self-similarity in the isotropic 
case leads to the Kolmogorov energy spectrum of distri- 
bution of energy among Fourier modes k, which reads 
E{k) ^ e^/'^fc"^/^, with e the energy injection rate. This 
law is rather well verifled in the atmosphere [5[ and in the 
solar wind [6|, |7[ , in laboratory [Sya] as well as when ex- 
amining numerical experiments [i3|- When considering 
the helicity spectrum H{k), where the helicity is the cor- 
relation between the velocity and its curl, the vorticity, 
one finds H{k) ^ k~^/^ (see for example |ll| - [l3j V How- 
ever, it was recently shown numerically and phenomeno- 



logically that, in the presence of rotation, other spectra 
emerge. This is the consequence of the effect of the iner- 
tial waves together with the fact that the helicity dom- 
inates the direct cascade to the small scales, the energy 
mostly undergoing an inverse cascade to the large scales 
[11, [H (see also [1^). 

Finding an underlying low-dimensional attractor to the 
complex behavior of turbulence is one way that such flows 
could be understood in simpler terms, as e.g., coherent 
structures with super-imposed noise. Another possibility 
is to resort to statistical mechanics and the dynamics of 
an ideal system with a finite number of modes. Equilib- 
rium statistical mechanics has proven valuable in that it 
could predict in many cases the direction of the transfers 
in wavenumber space, either to the small scales (direct) 
or to the large scales (inverse) p7H2l| . 

The Liouville theorem for fluid mechanics and mag- 
netohydrodynamics (when coupling the fluid equation to 
the temporal evolution of the magnetic induction) was 
flrst derived by T.D. Lee [13 (see also [23]). One can 
flnd a simple derivation of the Liouville theorem for two- 
dimensional (2D) turbulent flows in [21|. It expresses the 
incompressibility of the flow of dynamical variables in 
phase space (e.g., complex Fourier modes) using detailed 
conservation of quadratic invariants. Systems that can be 
described either by thermal equilibrium (zero-flux) or by 
non-equilibrium (finite fiux) dynamics can also be found 
outside the realm of fluid mechanics. An example is pro- 
vided by semi-conductor lasers described by the quantum 
Boltzmann equation |23| for which, as in many cases of 
optical turbulence, both a direct and an inverse cascade 
can be observed in the presence of pumping and dissipa- 
tion. Statistical equilibria have been derived for a variety 
of geophysical flows in the ideal case, in the presence of 
stratification [2^ or rotation [2^ [26 
and for quasi-geostrophic models 



for shear flows [27[ , 
2i, [2^. An equiva- 



TABLE I: Parameters of the runs: band of wavenumbers fco 
excited in the initial conditions, linear spatial resolution A*', 
rotation rate Q, initial relative helicity p, Rossby number Ro, 
and micro Rossby number Kui computed at f = (see text for 
definitions). At late times, R^^ is of order unity in all runs as 
a result of the thermalization. 
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lent problem in magnetohydrodynamics, when an exter- 
nal magnetic field is imposed, was studied in [30|. One 
of the concerns of these studies is the existence of quasi- 
invariants linked to the quasi bi-dimensionalization of the 
flow. In such case, the possibility of their condensation 
at the gravest mode available to the truncated system 
may indicate the existence of an inverse cascade in the 
forced-dissipative case (see e.g., |24|). 

However, the late-time statistical equilibrium of many 
of these systems seems to be independent of the geo- 
physical effects considered and only prescribed by the 
structure of the nonlinear terms in the Navier-Stokes 
equations. In addition, other invariants seem to play 
no significant roles. As an example, it was shown by 
Kraichnan [3l| that in the presence of helicity, the global 
ideal cquipartition equilibrium E(k) ~ k^ is only slightly 
altered at small scale by the presence of helicity, while 
the large-scales show no sign of condensation. In this 
context, it is worth asking what happens in the rotat- 
ing case, which can be viewed as a hybrid between the 
2D and the three-dimensional (3D) cases. In the forced- 
dissipative case, the small scales are fully 3D and isotropy 
is recovered at small scale, whereas the large scales are 
dominated by an inverse cascade of energy. The ideal 
truncated problem was tackled recently in (33] where it 
was shown that the establishment of the equilibrium is 
delayed in the presence of rotation in an anisotropic fash- 
ion. Furthermore, following analyses in [33|,|3J] and using 
a modal decomposition into wave, 2D, and 3D modes. 
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FIG. 1: Top: Energy E->,d (thick) and Eid (thin) (systemat- 
ically below Ezd) as a function of time in log- log coordinates 
for runs with initial conditions between fc = 2 and 4. Solid 
line: run Rl, dots: run R2, dash-dots: run Rll, all three 
with p = (see Table |I] and labels on curves giving the pair 
of values [fi, p]) and the linear resolution N . Finally, with 
p ~ 1 we have run R3 (dash line). Bottom: Ezd (thick) and 
E2D (thin) for runs with initial conditions between fc = 30 
and 40. All runs have p = and A'^ = 128. Solid line: run 
R14, dots: R17, dash: R19, dash-dots: R20 and dash-triple- 
dot: R21. The inset shows E2D for runs R17, R19, R20, and 
R21 with time in units of the inertial period rn (which goes 
from rn ~ 0.03 turnover times for run R17, to ~ 0.004 for 
run R21). 



[26| found that in non-helical flows in the limit of strong 
rotation (small Rossby number Ro), decoupling between 
inertial waves and 2D coherent structures obtains until a 
time t^, ^ Ro^ . After this time, coupling occurs but it 
can still be considered weak for long times and the system 
thus displays both 2D and 3D features. At this point, it 
is worth mentioning also that there is no consensus for 
the moment on whether perfect decoupling between the 
modes takes place for infinite rotation rate (see e.g., |3a|). 
We revisit here some of these results and extend our 
analysis to the helical case and to the dynamics of in- 
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FIG. 3: Temporal evolution of the enstrophy. Top: Runs with 
n = 16, initial conditions between fc = 2 and 4, for different 
spatial resolutions and values of the relative helicity. Bottom: 
Runs with initial conditions between fe = 30 and 40, N = 128, 
p = 0, and different values of fi. The inset shows the same 
curves with the time in logarithmic scale. 



FIG. 2: Top: E2D in simulations with initial conditions be- 
tween fc = 2 and 4, fixed spatial resolution (A^ = 128), and 
increasing fi. Non-helical runs R2 {Q. — 16), R4 (SI = 32), 
and R6 {0. = 64), and helical runs R3 {Q = 16), R5 {n = 32), 
and R7 (fi = 64) are shown. The arrows indicate the time 
when E2D decreases below 90% of its initial value, from left to 
right with increasing SI, above for helical runs and below for 
non-helical runs. The symbols are given in the figure, with 
[Tl, p] for /5 = in the bottom left, and for p = 1 in the 
top right. Bottom: E2D in simulations with initial conditions 
between k = 2 and 4, fixed rotation rate (Jl = 16), and in- 
creasing resolution (the higher, the thinner the line). Runs 
with zero helicity (dotted lines for runs R2, R8, and Rll), 
and with net helicity (dashed lines for runs R3, R9, and R12) 
are shown. 



termediate times during which an energy (and helicity) 
distribution can be observed reminiscent of dissipative 
dynamics, as shown in the non rotating case in [36l. |37|. 
In the next section, we describe briefly our procedure and 
move on in Sect. IIIII to discuss some of the predictions 
for early times in these systems, before thermalization is 
reached. Then, in Sect. IIVI we examine non-helical and 
helical flows initialized with perturbations either in the 
small or in the large scales. Section |V] presents our con- 
clusions. 



II. THE PROCEDURE 

We consider the equations for an incompressible and 
inviscid three-dimensional rotating fluid for the velocity 
fleld u, with V • u = 0. In the rotating frame, they read: 



ot 



(1) 



Here u; = V x u is the vorticity and V is the total pressure 
modifled by the centrifugal term, which is obtained self- 
consistently by taking the divergence of Eq. (HI), using the 
incompressibility condition, and solving for the resulting 
Poisson equation. We use a pseudo-spectral method in a 
{2tt)^ box with periodic boundary conditions. The code 
is fully parallelized using the message passing (MPI) li- 
brary j38l |39|; spatial resolutions in this paper are up 
to N^ = 512"^ with N points in each direction. The 
temporal scheme is a fourth order Runge-Kutta and the 
code uses the 2/3 dealiasing rule, resulting in a trun- 
cated set of modes from fcmin = 1 to a maximum wave 
number fcmax = N/3. The total energy E = (|up)/2 
and the total helicity H = (u-(£) are inviscid invari- 
ants preserved by the truncation [3ll |40| , with respective 
isotropic Fourier spectra E{k) and H{k); note that the 
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FIG. 4: Scaling of the time f** for a given scale to reach 
thermalization as a function of the imposed rotation; t^,^, ~ 
Q,'^''^ appears as a good fit. Different symbols correspond 
to different choices of the thermalization level (see text for 
details). Runs have initial conditions between k — 2 and 4, 
N = 128, and p = 0. 



Schwartz inequality implies \H{k)\/kE{k) < 1. The rel- 
ative hclicity will then be defined in the following as 



H{k)dk 




kE{k)dk j . 



(2) 



The rotation axis is taken to be in the z direction with 
fi = fJz, VL being the rotation frequency. Finally, the 
Rossby number is defined as usual as: 



Ro 
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2nLo 



(3) 



t/rms is the r.m.s. velocity taken to be of order unity, and 
Lq = 27r/fco is the characteristic scale of the initial condi- 
tions (fco in this expression corresponds to the minimum 
wavenumber with non-zero amplitude in the initial con- 
ditions). One can also define the micro Rossby number 
based on the r.m.s. vorticity 



Rou 



2n 



(4) 



measuring the intensity of the small-scale vorticity with 
respect to the imposed rotation. 

In order to examine long-time properties of the inviscid 
flow, computations were run for over 100 (in some cases 
up to 640) turnover times t = L/U-n-as in terms of a 
unit length L (in the following, time t in all figures is in 
units of r unless explicitly stated). The turnover time at 
the scale of the initial conditions is then given by tq ~ 
-^o/f^rms = {Lq/L)t and the inertial period is ro = Rotq. 

Two types of initial conditions are taken, one in the 
small scales, with the initial velocity field defined in the 
band of wavenumbers fco G [30,40], and one in the large 



scales, with fcg G [2,4]. The former case is to ensure 
that the effect of helicity is observable in the absolute 
equilibrium derived by Kraichnan [3l| , whereas the latter 
is to examine a sufficiently resolved energy spectrum in 
the direct transfer to small scales. Also, the two sets of 
initial conditions allow us to study the dependence with 
scale of the time to reach the statistical equilibrium. The 
initial amount of helicity is controlled using the method 
given in [4l|, coupling two random fields by imposing a 
prescribed angle (in physical space) between them. Some 
of the properties of the runs described in this paper are 
summarized in Table ID 

Assuming a truncated system of modes and consid- 
ering the two quadratic invariants, the energy E = 
J E{k)dk and the helicity H ~ J H{k)dk, the macro- 
canonical probability density function is proportional to 
Cexp [~aE ~_BH]. This results in the ideal solutions for 
the spectra j31| : 



E{k)^ 



2afc2 
D{k) 



H{k) = 



2/3fc4 



(5) 



with 



D{k) ^ a^ - P^k^ 



In these solutions, a and /3 can be seen as temperatures 
associated with the two invariants E and H . constrained 
by a > and |/3|fcmaa; < ol (for given initial energy and 
helicity, these equations can be solved for a and /3; see 
[37|). Note that helicity and thus /3 are not positive def- 
inite: they are pseudo-scalars, changing sign when one 
goes from a right-handed to a left-handed frame of ref- 
erence. Since H{k)/kE{k) ~ fc in Eq. ([5]), eddies at a 
smaller scale are more helical, contrary to the dissipative 
case for which this ratio is observed to vary as 1/fc in the 
non rotating case, thereby indicating a slow ^ k~^ re- 
covery of mirror-symmetry as fc — >■ oo. This discrepancy 
in the relative behavior of spectra may be the origin of 
persistence of helical motions in the dissipative range of 
viscous fiows [H, 113 ■ With /3 = when iJ = 0, one 
recovers in the ideal case the equipartition of energy be- 
tween Fourier modes, with a number of modes per spher- 
ical shell of radius fc and unit width proportional to fc^ . 
In the presence of solid-body rotation the Coriolis force, 
being conservative and linear in the velocity, does not al- 
ter the long-time statistical equilibria given in equation 
([S]), and thus the isotropic distribution of energy is also 
unaltered. 

The numerical conservation of energy and helicity is 
therefore essential when performing ideal truncated runs 
and thus we have systematically monitored their evolu- 
tions. For example, in runs with Vt ~ 128, we had to 
resort to a time-step 32 times smaller than in runs with 
r2 = (resulting, e.g., in Ai = 5 x 10"'^ for the former 
and At « 1.5 x lO"'' for the latter when TV = 128). In 
such a case, the energy is conserved to within 0.5% after 
t w 100 for il = 128, and better still for smaller rotation 
rates at the same resolution (e.g., at 0.02% for il ~ 6A at 
the same time). 
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FIG. 5: Time evolution of the isotropic energy spectrum. 
Above are runs with A'^ = 128: (a) Rl with Jl = and p = 0, 
and (b) R2 with fi = 16 and p = 0; in (a) and (b) times are 
t = (solid), t = 4 (dotted), f = 10 (dashed), t = 40 (dash- 
dotted), i = 80 (dash-triple-dotted), t = 150 (long dashed), 
and t — 320 (dotted). Below are runs with N — 256: (c) 
Rll with n = 16 and p = 0, and (d) R12 with fl = 16, and 
p « 1; in (c) and (d) times are (with the same line style order 
as above) f = 0, f = 4, f = 8, i = 16, i = 40, f = 80, and 
t — 160. Different slopes are indicated only as references. All 
runs build a k^ tail over time. 
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FIG. 6: Energy flux n(fc) in runs Rl (a) and R2 (b). Symbols 
for different times are as in Fig. [Sfa-b). 



III. IS THERE A LARGE-SCALE IDEAL 
CONDENSATE IN THE ROTATING CASE? 

As mentioned previously, an inverse energy cascade has 
been observed in several numerical simulations of driven 
dissipative rotating turbulence whereas the (isotropic) 
statistical equilibria show no sign of large-scale condensa- 
tion. We will show that in our simulations, the obtained 
equilibria are independent of the resolution used. How- 
ever, very long integrations are needed in the presence 
of rotation and helicity to reach these final thermalized 



states. As a result, we next resort to a series of direct 
numerical simulations (DNS) at moderate resolution and 
computing for long times to examine further this prob- 
lem. 

It is useful to decornpose the total velocity field in the 
following manner [SJ, |42, |4j] . We define a preferred di- 
rection z, parallel to the rotation axis. A general three- 
dimensional fiow u can then be written as the sum of 
a strictly two-dimensional velocity u^(A:^,A:|| = 0) of 
energy E±, a vertical component Uz{k±,ku — 0) of en- 
ergy Ez , and the remaining three-dimensional component 
U3(/c±,fc|| 7^ 0) of energy E^y). The sub-indices _L and || 
refer to the rotation axis. The total kinetic energy is: 



E ^ ^ f \u{x)\^(fx= I E{k)dk 



= E: 



3D 



■E^ 



Ez — E; 



3D 



E- 



2D, 



(6) 



where E2B 



E± + Ez is the total energy in the modes 
with fc|| = 0. Note that there is a one to one corre- 
spondence [431 between E^^ and e(fc_L,fc|| ^ 0), E±_ oc 
[e — ZpA](fc±,fc|| = 0), and Ez oc [e + 2'pA](fcx,fc|| = 
0), where e(fc_L,fc||) is the two-dimensional axisymmet- 
ric energy spectrum, and ZpA{k_\_^k\\) is the polarization 
anisotropy. The spectra can be expressed in Cartesian co- 
ordinates (fcj_,fc||), or in polar coordinates (fc,cos0) where 
9 is the angle between the vector k and the vertical axis. 
These scalar spectra e and ^pa, together with the helic- 
ity spectrum, are the dynamical variables used in studies 
of rotating turbulence using spectral closures [ijl . 

Bourouiba [261 argues, by considering only resonant 
triads, that for small enough Rossby number, the 2D and 
3D modes are effectively decoupled for times i < t, . This 
implies separate conservation properties of £'30 , E±_ , and 
Ez- In that case, a different statistical equilibrium to the 
one described in Sect. |ll] can be derived, which leads to 
separate predictions for the horizontal dynamics of the 
ffow. For late times (i > t^), the isotropic statistical 
equilibrium is recovered. 

A non-ideal example where decoupling between 2D and 
the 3D modes can be observed may be given by freely de- 
caying rotating flows. In [4^, it was found that in this 
case the 2D and 3D energies decay at different (and in- 
dependent) rates. If we follow the estimation of coupling 
being effective after the time i* ~ Ro~ [431, ^^'^ taking 
into account that the Rossby number itself decays as t^/^ 
according to the study in [45|, we conclude that i* ~ t. 
In other words, t* may never be reached. It is not clear, 
however, that in the forced case where the Rossby num- 
ber is maintained at a constant value, such a decoupling 
can remain permanently [46|. 



IV. NUMERICAL RESULTS 

A. Temporal dynamics properties 

We first consider an initial condition concentrated in 
the large scales (runs Rl to R13). An isotropic excitation 
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FIG. 7: Top: Isotropic energy spectrum at f = 5 (at the be- 
ginning of the intermediate phase when thermalization begins 
at small scale) in runs without helicity and with Q. = 16: run 
R2 (solid), Rll (dotted), and R13 (dashed), with increasing 
resolution from N = 128 to A'^ = 512. Bottom: energy spec- 
trum compensated by k~ at t — 10 in 256 runs RIO (solid 
line, SI = 8 and p = 0), Rll (dotted, fl = 16 and p = 0), and 
R12 (dashed, n = 16 and p = 1). 



with random phases and fast decaying energy spectrum 
is put at t = between k = 2 and 4 in Fourier space, such 
that Urms ~ 1 and tq ~ 3.1. The temporal evolution of 
the energies, specifically E^d and E2D (see Eq. ([6])) are 
displayed in Fig. [T] for rotation rates il = and il = 16, 
and for different spatial resolutions ranging from N = 
128 to 256. Simulations with and without helicity are 
considered. 

In the 128"^ run with il = and p = 0, E^i^ grows 
fast to « 0.49 and E2B decays to « 0.01 after t « 10. 
These values correspond to the isotropic statistical equi- 
librium, where all the modes have thermalized and have 
the same energy per mode (-E3D is larger as there are 
more 3D modes than modes in the fcy = plane). In the 
128'^ run with fl = 16, these thermalized values are only 
reached after t « 300. In the presence of rotation not 
only the decay towards the isotropic statistical equilib- 
rium is slowed down; at early times (up to t « 4 « 1.3ro) 
a transient regime develops for which E2D and E^u seem 
to be independently conserved. Details of this regime 
for different Rossby numbers can be found in [26|. Here 
we consider the effect of resolution and helicity in the 
time when the quasi-conservation breaks down. To first 



FIG. 8: Two dimensional energy spectrum e{k, cos 6) in run 
R13 at t = 10, for e = (dashed), 6 = 7r/4 (dotted), and 
9 = 7r/2 (solid). A slope of —2 is shown as a reference. 



illustrate this, a 256'^ run with O = 16 is also shown in 
Fig. [TJ The change in resolution does not seem to affect 
the quasi-conservation of ,£20 and E^d , although the fol- 
lowing decay towards the isotropic statistical equilibrium 
takes place at a slower pace at larger resolution. A simi- 
lar result was obtained for a 512'^ simulation (not shown). 
Finally, the presence of helicity in a 128'^ run with il = 16 
increases (by close to a factor 2) the time during which 
E2D and E3D can be considered as quasi-invariants. 

The effect of resolution and helicity in the quasi- 
conservation can be further explored considering runs R2 
to R12. Runs R2 to R7 span fi from 16 to 64 with initial 
relative helicity p of or 1 at fixed resolution (A^ = 128). 
The time evolution of E2D in these runs is shown in 
Fig. [5J As ri is increased, the time during which E2D 
(and E^o) can be considered quasi- invariant increases, 
and the differences between the non-helical and helical 
cases become more evident. In all cases, helicity extends 
the time of quasi-conservation. At the end of the quasi- 
conservation phase, E2D in all runs decreases with nearly 
the same rate independently of fi, and depending only on 
whether the runs have helicity or not. 

Figure [2] also shows the effect of changing resolution. 
Runs with fixed Q are shown with and without helic- 
ity, while changing TV from 64 to 256. The time during 
which E2D is quasi-invariant seems insensitive to the lin- 
ear resolution, while the rate of change after the quasi- 
conservation (as well as the final value reached at late 
times) depends on N . As a result, high resolution is not 
necessary to study the early time (or, as shown next, the 
very late time) dynamics of the ideal truncated systems. 
Simulations at higher resolution in Table |T] will mostly be 
useful to consider the spectrum that develops at interme- 
diate scales during the transition from the quasi-invariant 
regime to the isotropic thermalized regime. 

The case of initial conditions concentrated in the small 
scales (runs R14 to R21) is also shown in Fig. [TJ Again 
an isotropic excitation at t = is considered, but now 
between wavenumbers A: = 30 and 40, with Urms = 1 
and To ~ 0.2. Except for very early times, the E2D and 



£^30 energies remain constant and close to « 0.49 and 
« 0.01 respectively. However, unlike the previous runs 
and the cases considered in [2a] , the constancy of E2B and 
i^sD should not be interpreted as a conservation of quasi- 
invariants (note it also takes place for f2 = 0). With the 
initial conditions at small scales, the system goes faster to 
the isotropic thcrmalized equilibrium (the turnover time 
is shorter in this case), and the constancy of E2B and 
E3D is simply the result of all modes having thcrmalized 
with the same energy. Once the thermal equilibrium is 
achieved, the solution is flux-less and thus no exchange 
of energy takes place. An approximate conservation of 
the quasi-invariants may happen at much earlier times, 
as shown in the inset of Fig. [TJ Runs R20 with fi = 64 
and R21 with H. = 128 show an early increase of E21) 
followed by approximate constancy, to later decay to the 
isotropic thermalized value. 

As a result, the time to reach equilibrium is delayed in 
the presence of either rotation or hclicity, but the evolu- 
tion of such flows drives them to the three-dimensional 
equilibrium given in Eq. ([S|) as will be shown in more 
detail later. This delay can also be observed in the time 
evolution of the enstrophy. In Fig. [3] we contrast the 
effects on its temporal evolution of the scale of the ini- 
tial conditions, of the initial relative helicity, and of the 
spatial resolution. In the runs with initial conditions 
between k = 2 and 4, (w^) grows fast and reaches its 
saturated value (which corresponds to the isotropic ther- 
malized state) after t w 100. Increasing the resolution 
does not affect this time, although it substantially in- 
creases the saturation value of the mean square vorticity 
(as more wavenumbers are available for the thermaliza- 
tion). On the other hand, the presence of helicity slows 
down the decay towards the thermalized state. Similar 
results are obtained in the runs with initial conditions 
between fc = 30 and 40. In this case, (w^) decreases 
with time as the initial enstrophy, concentrated at small 
scales, is distributed among all wavenumbers as the sys- 
tem thermalizes. These runs also illustrate how the time 
to reach the thermalization increases as H. is increased. 

Defining an intermediate wavenumber between fcmin 
and fcmax, we can compute the time i*, at which ther- 
malization has reached that wavenumber, in the sense 
that its amplitude does not evolve except for thermal 
fluctuations around the mean after that time. In prac- 
tice, this is equivalent to measuring the time at which a 
certain value of (i^^) is reached in Fig. |31 Fig. U displays 
the scaling of that time as a function of imposed rotation 
rate; a J7^/"' range appears clearly. We have checked that 
this scaling is insensitive to the choice of scale at which 
it is computed. The different symbols in the figure cor- 
respond to different choices of thermalization level, with 
diamonds corresponding to (w^) within 1% of the ther- 
malized value, triangles to 2%, squares to 3%, and crosses 
to 4%. Note this time should not be confused with the 
time t* for coupling of 2D and 3D modes that results in 
the end of the conservation of quasi-invariants. 
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FIG. 9: (a) Isotropic energy spectrum, (b) reduced perpen- 
dicular spectrum, and (c) reduced parallel spectrum at f = 44 
in non-helical runs R14 (solid, ^ = 0), R17 (dotted, fi = 16), 
R19 (dashed, fi = 32), R20 (dash-dotted, fi = 64), and R21 
(dash-triple-dotted, Q, — 128). Slopes are given as references. 



B. Large-scale dynamics and early time evolution 

As it will be confirmed in the next subsection, the ideal 
flow ends up in a state of cquipartition of energy as pre- 
dicted by T.D. Lee [l7| (or its generalization to helical 
flows [31|). However, the path to get there is interest- 
ing as it uncovers a "dissipative" transient that displays, 
in the non-rotating case, a Kolmogorov spectrum. Al- 
though the total energy is conserved, the large wavenum- 
ber k^ tail - established at intermediate times as a partial 



thermalization - plays the role of a turbulent viscosity 
for the large scales, which can undergo a standard turbu- 
lence dynamics. This was shown both for the non-helical 
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FIG. 10: (a) Isotropic spectra, (b) reduced perpendicular 
spectra, and (c) reduced parallel spectra of the energy (thin) 
and of the helicity (thick) at f = 55 in helical runs R15 (solid, 
SI = 0) and R18 (dotted, Q. = 16). Slopes are given as ref- 
erences. Helicity spectra are normalized by the total helicity. 
In (a) , the isotropic prediction of Kraichnan [31[ is shown in 
dashed lines, with values of a and /3 obtained from solving 
Eq. dSJ. 



case [36[ and for the helical case [37[. We thus exam- 
ine this temporal behavior in the presence of rotation, 
by considering the runs with initial conditions at large 
scales in the band of wavenumbcrs between fc = 2 and 
4. Note that this corresponds to the transition of the 
system from the quasi-invariant regime to the isotropic 
statistical equilibrium. 

The time evolution of the isotropic energy spectrum is 
shown in Fig. [5] for a non-helical non-rotating run (run 
Rl), two rotating non-helical runs (R2 and Rll, at dif- 



ferent resolutions), and a rotating helical run (R12). Sev- 
eral times are displayed in order to show the progressive 
advance of thermalization towards large scales and the 
dynamics at large scale. Run Rl at early times devel- 
ops a spectrum compatible with a ^ fc~^'^ law, and as 
the thermalized range widens the spectrum goes progres- 
sively to the ~ fc^ solution. For non-rotating flows, this 
behavior was studied before (we refer the reader for de- 
tails to Refs. [H,!!! 113, 113). In the rotating case (Run 
R2), a dissipative-like transient also develops as the sys- 
tem transits from the regime with quasi-conserved E2T) 
and £^3D to the isotropic thermalized equilibrium, but 
the large-scale spectrum is steeper in this case. As will 
be shown next, the steeper spectrum seems to be con- 
sistent with a ~ fc~^ law. As in the non-rotating case, 
the small scales evolve towards the isotropic ^ fc^ solu- 
tion, but the time for the large-scales to reach thermaliza- 
tion is substantially increased in the presence of rotation. 
Changing spatial resolution does not seem to affect these 
results, but the presence of helicity makes the large-scale 
spectrum even steeper (closer to ^ fc~^. compare (c) and 
(d)). 

Figure IH gives the energy fluxes at different times for 
the 128'^ non-helical runs with J7 = and 16 (runs Rl and 
R2). These fluxes evolve towards zero as time elapses, 
corresponding to zero-flux statistical equilibria solutions, 
but at intermediate times they show clear energy ex- 
changes between Fourier modes. In the case with i7 = 0, 
a range with positive flux (which can be associated to 
the ^ k~^/^ range) is observed before the flux goes to 
zero at all scales. In the presence of rotation, besides a 
range with approximately constant direct flux, a change 
of sign of the flux at large scales is observed which indi- 
cates a tendency toward a large-scale condensation, as- 
sociated to the inverse cascade of energy in the forced- 
dissipative case. However, the final state lacks such con- 
densate as already stated, and the system evolves towards 
zero-flux equilibria (although at later times than in the 
non-rotating case). 

The effect of resolution and of helicity in the large-scale 
energy spectrum is further illustrated in Fig. [71 which 
shows the spectrum at i = 5 in runs with $7 = 16 and 
p = for increasing spatial resolution up to A^ = 512. At 
this early time little thermalization '^ k"^ is present yet. 
All runs show spectra compatible with a k~^ law at large 
scales although a knee followed by a slightly shallower 
spectrum is observed in the higher-resolved runs. As will 
be shown next, this knee in the spectrum is the result of 
anisotropics developing in the flow. Figure [7] also gives 
the energy spectrum compensated by k~^ at f = 10 in 
runs with 256'^ grid points and for various rates of rota- 
tion and of relative helicity. The change of Vt does not 
affect the power law, although in the presence of helicity 
a steeper energy spectrum obtains, in agreement with the 
results shown in Fig. [5j Also, a slightly shallower spec- 
trum than fc^ obtains in the smallest scales, indicating 
again a slowing-down of the dynamics. 

Both in the non-helical and helical cases in the presence 
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FIG. 11: Time evolution of the energy in the 2D mode k — 
(1, 1, 0) after perturbing the thermalized value by multiplying 
its amplitude by 10 in runs R14 (thin gray, Q = 0) and R17 
(thick black, Q — 16). The inset shows the time evolution for 
late times; note the oscillations in run R17, corresponding to 
Rossby waves. 



of rotation, a spectrum different from Kolmogorov is ob- 
tained, with an inertial index between —3 and —4. Such 
spectra also differ from the ones found in driven dissipa- 
tive studies of rotating turbulence, namely E{k) ~ fc~^ 
or a bit steeper in the helical case. For non- helical rotat- 
ing turbulence, a k~^ spectrum has been advocated using 
a spectral closure [35| because of the distribution of en- 
ergy for different angles between the wavevector k and 
the rotation axis, the classical fc~^ spectrum being found 
only for modes with ku « 0. Indeed the present simu- 
lations confirm this prediction. The axisymmetric two- 
dimensional energy spectrum e{k,cos9) [J, |49| is shown 
in Fig. [5] for different values of 6 for run R13 at i = 5. 
The differences in the spectra for different directions in- 
dicate the development of anisotropics, with most of the 
energy in modes with fcii = {9 = 0) and following 
e{k,cos9 = 1) '-^ k~^. As 9 increases towards 6 — tt/2 
(corresponding to wave vectors k = fc||z) the energy con- 
tent decreases and the spectrum becomes steeper. Note 
that an isotropic fc~^ spectrum can be recovered using a 
phenomenological analysis taking into account the slow- 
ing down of nonlinear transfer due to inertial waves, still 
assuming isotropy [50, IMl . The extension of this argu- 
ment to the anisotropic case is straightforward and leads 



to E{k±,k\i) 



k-'/'k; 



1/2 



However, the anisotropic 



spectra observed here in the transient phase are steeper 
and are more likely attributed to an angular dependence 
as advocated in 13511 . 



C. Isotropic statistical equilibrium 

At late times, all the runs (see e.g.. Fig. [5]) reach the 
statistical equilibrium with zero flux and energy spec- 
trum ~ k^ with no sign of condensation at large scales. 
To study this later regime in more detail, we focus now 
on the runs with initial conditions at small scales (runs 
R14 to R21), as these runs reach thermalization faster. 
Similar results are obtained at late times in runs Rl to 
R13. 

We show in Fig. ^ the energy spectra for different ro- 
tation rates (see caption) at t ^ 44 ^ 220to. Note the 
evolution of all runs toward a pure fc^ spectrum, as ex- 
pected from the isotropic equations ([5]) . Only the modes 
at the smallest wavenumbers depart from the equilib- 
rium, as the time it takes for these modes to thermal- 
ize increases as fl increases. This state is isotropic, as 
can be confirmed, e.g., by studying the "reduced" par- 
allel and perpendicular energy spectra E(ku) and E{k±) 
(see e.g., [5j, (U). The spectrum E{k\\) is obtained by 
integrating the energy over all modes with same wave 
vector fell (i.e., integrating over planes in spectral space), 
and E{k±) is obtained by integrating the energy over 
all modes with same wave vectors fcj_ = (fc^ -I- fc^)^^^ 
(i.e, integrating in cylindrical shells). The perpendicu- 
lar spectrum shows a ^ fc dependence (indicated by the 
solid line for convenience) , as expected for equilibria with 
integration over wavenumbers in one direction only (see 
also [201). The parallel spectrum is flat, again consis- 
tent with equipartition of energy between all modes, as 
in that case all planes in Fourier space contain the same 
amount of energy. The decrease of both reduced spec- 
tra at large wavenumber is due to the fact that, in the 
anisotropic case, integrations are performed in planes or 
cylinders but undergo the spherical truncation imposed 
in the computation (associated to the dealiasing). This 
results in a deficit of modes at large |k| for anisotropic 
binning of Fourier space. 

When helicity is present, the absolute equilibria en- 
semble are modified as found in |3l| . with a behavior 
close to k^ at low wavenumber but departing from it at 
high k. The result is illustrated in Fig. [TUl which shows 
the reduced energy and helicity spectra in runs R15 and 
R18 at i = 55 = 275to. The reduced helicity spectra are 
constructed following the same procedure as for the en- 
ergy. Independently of the rotation rate, the spectra are 
in good agreement with the isotropic theoretical predic- 
tions given in [31| . The delay for reaching the asymptotic 
state is slightly larger in the helical case at the same ro- 
tation rate, as can be expected from the fact that for a 
fully helical Beltrami flow, the Lamb vector u x a; w 
until the initial flow develops. 

In light of these results, it may seem that the final 
equilibrium is exactly the same for all the runs, inde- 
pendently of the presence of rotation. However, there 
are differences in the coupling between modes in the ro- 
tating case even after thermalization has been achieved. 
To show this, the final states of runs R14 and R17 were 
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perturbed by multiplying by 10 the amplitudes of the 
modes k = (1, 1,0) (a 2D mode) and k = (10, 10, 10) (a 
3D mode). Return to thermalization of these modes took 
different times depending on the value of fi, more clearly 
in the case of the 2D mode (see Fig. fTTj) . Moreover, once 
the thermalized amplitude was reached, the modes in the 
run without rotation fluctuated in time as white noise, 
while the same modes in the run with rotation showed 
clear oscillations with the frequency of the Rossby waves. 



V. DISCUSSION AND CONCLUSION 

We have found in the case of helical rotating flows 
a counter-example to the otherwise well-tested conjec- 
ture that statistical ensemble equilibria are good predic- 
tors for the small-scale and large-scale behavior of tur- 
bulent flows. In three-dimensional helical rotating flows, 
the equilibrium does not have a condensate with a low 
wavenumber accumulation of excitation, whereas the dis- 
sipative dynamics show a clear inverse cascade. More- 
over, the spectra computed for long times for inviscid 
truncated systems do not show any anisotropization of 
the large scale. The limitation of the conjecture that sta- 
tistical ensemble equilibria are good predictors for sniall- 
and large-scale behavior was already identified for the 
case of non-helical rotating flows in [2y, |32| . 

One is led to speculate as to what is the origin of this 
discrepancy; we attribute it to the way the nonlinearity of 
the system is modified in the presence of inertial waves, as 
studied from numerical data, e.g., in [Tj, [5^. In rotating 
turbulence, the development of anisotropy and inverse 
transfer of energy results from the resonant and near- 
resonant interactions between three-dimensional modes 
and modes in the slow manifold with fcii = 0, per virtue of 
the linear Coriolis term [42], |4J| . Energy in three dimen- 
sional modes is transferred by a subset of the resonant 
interactions to modes with smaller vertical wavenumber. 
It is therefore to be expected that the inverse transfer of 
energy in a rotating flow will be different from the two- 
dimensional case (without rotation) insofar as resonances 
play an essential role. However, the Liouville theorem is 
not modified by the linear terms, and therefore the statis- 
tical equilibrium remains isotropic and without conden- 
sates after rotation is introduced, except for early times 
when the energy in the 2D modes is quasi-conserved. It 
remains puzzling however to see that the waves, stem- 
ming from a linear term in the dynamical equation, per- 
turb a nonlinear effect when those end up dominating 
the dynamics for sufficiently long times. This paper at- 
tempted at elucidating this behavior examining the role 
that the second invariant, the helicity, has on inviscid 
dynamics since it is known that it affects the driven dis- 
sipative case, making the spectra different 11 5| . 1551 from 
what is expected in the non- helical case 50|, |51| . 

It has been argued [2^, [2y| that, since in the case 
of very strong rotation, the flow evolves toward a two- 
dimensional state in which two-dimensional quadratic in- 



variants as the enstrophy are recovered, this can justify 
the existence of an inverse cascade in the driven dissi- 
pative case. Indeed, we can confirm the results of [2g| 
showing quasi-conservation of two-dimensional invariants 
for early times. The introduction of helicity further in- 
creases the time in which this approximation holds, but 
keeps as the equilibrium state that of the non-rotating 
case. However, independently of the presence of helic- 
ity. once the isotropic thermalized equilibrium has been 
reached, differences in the evolution of the modes can be 
identified. While in the absence of rotation the amplitude 
of the modes evolve as white noise, traces of rotation and 
the associated coupling between modes can be identified 
in the rotating case. Perturbation of any mode requires 
long times to recover thermalization in the latter case, 
and time oscillations associated to waves are observed in 
the amplitude of the modes. 

The results are insensitive to the spatial resolution 
used, at least in the range of values explored here. The 
time during which E2U and £"30 are independently con- 
served seems insensitive to the resolution. Properties of 
the late time isotropic equilibrium are also unchanged 
as resolution is increased. However, resolution allows us 
to resolve better the spectrum that develops during the 
transient. Our results show, in a gree ment with previ- 
ous claims for non-rotating flows [3y|, that during the 
transition from the approximate conservation of quasi- 
invariants to the thermalized regime, a dissipative-like 
spectrum arises. There is a dispute as to what the energy 
spectrum for rotating flows may be. even when putting 
aside the effect helicity may have [l 5 | .[55l . Weak turbu- 
lence theory predicts a k ~^ sp ectrum |56| . consistent with 
dimensional analysis [50, ^M (note that this theory finds 
at lowest order a decoupling of the 2D and 3D modes). 



But it is claimed in 



that such a spectrum applies 
= and that, as one takes into 



only in the vicinity of ku 
account the contribution from all parallel wavenumbers, 
the resulting isotropic spectrum is oc k~^. Note that this 
power-law should not be viewed as the classical Kraich- 
nan spectrum for two-dimensional Navier-Stokes [2l| , but 
rather as being the isotropic result of an anisotropic dis- 
tribution of modes taking into account the fcy dependence 
of the dynamics. 

Essential to the argument is that the decoupling of 2D 
and 3D motions on the basis of resonances [i^l is not com- 
plete because of the differential energy between purely 2D 
and purely 3D vertical velocities in the fc|| case (the polar- 
ization anisotropy Zpa)- The results presented here seem 
to favor this spectrum, with an isotropic scaling '^ k~^ 
which obtains from an anisotropic two-dimensional spec- 
trum e(fc,cos6') which scales for cos 6* = 1 as fc~^ (in 
the presence of helicity, the isotropic spectrum is even 
steeper). As a result, the ideal runs seem therefore to be 
a useful tool to study the high-rotation case in spite of 
the limitations of the statistical ensemble predictions. It 
is unclear for the moment why in viscous simulations of 
non-helical rotating flows the isotropic spectrum scales 
as k~^. It may be the case that dissipation plays an es- 
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scntial role at sufficiently small scale, an effect already 
noted in [53] in the decaying case. 
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